Experimental and In Silico Analysis of TEM ß-Lactamase Adaptive Evolution.

Multiple mutations often have non-additive (epistatic) phenotypic effects. Epistasis is of fundamental biological relevance but is not well understood mechanistically. Adaptive evolution, i.e., the evolution of new biochemical activities, is rich in epistatic interactions. To better understand the principles underlying epistasis during genetic adaptation, we studied the evolution of TEM-1 ß-lactamase variants exhibiting cefotaxime resistance. We report the collection of a library of 487 observed evolutionary trajectories for TEM-1 and determine the epistasis status based on cefotaxime resistance phenotype for 206 combinations of 2-3 TEM-1 mutations involving 17 positions under adaptive selective pressure. Gain-of-function (GOF) mutations are gatekeepers for adaptation. To see if GOF phenotypes can be inferred based solely on sequence data, we calculated the enrichment of GOF mutations in the different categories of epistatic pairs. Our results suggest that this is possible because GOF mutations are particularly enriched in sign and reciprocal sign epistasis, which leave a major imprint on the sequence space accessible to evolution. We also used FoldX to explore the relationship between thermodynamic stability and epistasis. We found that mutations in observed evolutionary trajectories tend to destabilize the folded structure of the protein, albeit their cumulative effects are consistently below the protein's free energy of folding. The destabilizing effect is stronger for epistatic pairs, suggesting that modest or local alterations in folding stability can modulate catalysis. Finally, we report a significant relationship between epistasis and the degree to which two protein positions are structurally and dynamically coupled, even in the absence of ligand.

Gemcitabine and cisplatin plus durvalumab with or without tremelimumab in chemotherapy-naive patients with advanced biliary tract cancer: an open-label, single-centre, phase 2 study.

BACKGROUND: Immunotherapies have shown clinical activity in patients with advanced biliary tract cancer, for which outcomes remain poor despite standard of care treatment with gemcitabine and cisplatin. We aimed to evaluate gemcitabine and cisplatin plus durvalumab with or without tremelimumab as first-line treatment in patients with advanced biliary tract cancer. METHODS: This open-label, single-centre, phase 2 study was conducted at Seoul National University Hospital. Eligible patients were treatment-naïve adults aged 18 years or older with histologically proven unresectable or recurrent biliary tract cancer, at least one measurable lesion based on the Response Evaluation Criteria in Solid Tumors (version 1.1), an Eastern Cooperative Oncology Group performance status of 0 or 1, life expectancy of 12 weeks or longer, and adequate healthy organ and bone marrow function. Initially, all patients received one 3-week cycle of gemcitabine (1000 mg/m2) and cisplatin (25 mg/m2) on day 1 and 8 followed by gemcitabine and cisplatin plus durvalumab (1120 mg) and tremelimumab (75 mg) on day 1 of each cycle, starting with the second cycle (chemotherapy followed by chemotherapy plus durvalumab and tremelimumab group). Following protocol amendment, patients were recruited to receive gemcitabine and cisplatin plus durvalumab, starting on day 1 of the first cycle (chemotherapy plus durvalumab group) or gemcitabine and cisplatin plus durvalumab and tremelimumab also from day 1 of the first cycle (chemotherapy plus durvalumab and tremelimumab group) in parallel and allocated using a random block method. Assessors and patients were not masked to the treatment group. The primary endpoint was objective response rate, assessed in the efficacy population (ie, patients who were treated at least until the first tumour response assessment). This study is registered with ClinicalTrials.gov, NCT03046862 (active). FINDINGS: Between March 2, 2017, and Feb 13, 2020, 128 patients were enrolled (32 in the chemotherapy followed by chemotherapy plus durvalumab and tremelimumab group, 49 in the chemotherapy plus durvalumab group, and 47 in the chemotherapy plus durvalumab and tremelimumab group). Four patients (two in the chemotherapy followed by chemotherapy plus durvalumab and tremelimumab group and two in the chemotherapy plus durvalumab group) were excluded and 124 were evaluable for tumour response. The median duration of follow-up was 48·2 months (IQR 41·5-49·4) for the chemotherapy followed by chemotherapy plus durvalumab and tremelimumab group, 26·6 months (19·0-27·9) for the chemotherapy plus durvalumab group, and 24·2 months (20·7-31·7) for the chemotherapy plus durvalumab and tremelimumab group. 82 (66%) of 124 patients achieved an objective response (15 [50%] of 30 in the chemotherapy followed by chemotherapy plus durvalumab and tremelimumab group, 34 [72%] of 47 in the chemotherapy plus durvalumab group, and 33 [70%] of 47 in the chemotherapy plus durvalumab and tremelimumab group). The most common grade 3 and 4 adverse events were decreased neutrophil count (67 [53%] of 126), anaemia (50 [40%]), and decreased platelet count (24 [19%]), with no unexpected safety events. No adverse events leading to discontinuation or death occurred. INTERPRETATION: Gemcitabine and cisplatin plus immunotherapy showed promising efficacy and acceptable safety in patients with biliary tract cancer. Gemcitabine and cisplatin plus durvalumab are being evaluated in the phase 3, TOPAZ-1 study (NCT03875235) as first-line treatment in patients with advanced biliary tract cancer. FUNDING: AstraZeneca; National Research Foundation of Korea (Grant No. 2021R1A2C2007430).

Validation of a Circulating Tumor DNA-Based Next-Generation Sequencing Assay in a Cohort of Patients with Solid tumors: A Proposed Solution for Decentralized Plasma Testing.

BACKGROUND: Characterization of circulating tumor DNA (ctDNA) has been integrated into clinical practice. Although labs have standardized validation procedures to develop single locus tests, the efficacy of on-site plasma-based next-generation sequencing (NGS) assays still needs to be proved. MATERIALS AND METHODS: In this retrospective study, we profiled DNA from matched tissue and plasma samples from 75 patients with cancer. We applied an NGS test that detects clinically relevant alterations in 33 genes and microsatellite instability (MSI) to analyze plasma cell-free DNA (cfDNA). RESULTS: The concordance between alterations detected in both tissue and plasma samples was higher in patients with metastatic disease. The NGS test detected 77% of sequence alterations, amplifications, and fusions that were found in metastatic samples compared with 45% of those alterations found in the primary tumor samples (p = .00005). There was 87% agreement on MSI status between the NGS test and tumor tissue results. In three patients, MSI-high ctDNA correlated with response to immunotherapy. In addition, the NGS test revealed an FGFR2 amplification that was not detected in tumor tissue from a patient with metastatic gastric cancer, emphasizing the importance of profiling plasma samples in patients with advanced cancer. CONCLUSION: Our validation experience of a plasma-based NGS assay advances current knowledge about translating cfDNA testing into clinical practice and supports the application of plasma assays in the management of oncology patients with metastatic disease. With an in-house method that minimizes the need for invasive procedures, on-site cfDNA testing supplements tissue biopsy to guide precision therapy and is entitled to become a routine practice. IMPLICATIONS FOR PRACTICE: This study proposes a solution for decentralized liquid biopsy testing based on validation of a next-generation sequencing (NGS) test that detects four classes of genomic alterations in blood: sequence mutations (single nucleotide substitutions or insertions and deletions), fusions, amplifications, and microsatellite instability (MSI). Although there are reference labs that perform single-site comprehensive liquid biopsy testing, the targeted assay this study validated can be established locally in any lab with capacity to offer clinical molecular pathology assays. To the authors' knowledge, this is the first report that validates evaluating an on-site plasma-based NGS test that detects the MSI status along with common sequence alterations encountered in solid tumors.

Multiregion whole-exome sequencing of intraductal papillary mucinous neoplasms reveals frequent somatic KLF4 mutations predominantly in low-grade regions.

OBJECTIVE: Intraductal papillary mucinous neoplasms (IPMNs) are non-invasive precursor lesions that can progress to invasive pancreatic cancer and are classified as low-grade or high-grade based on the morphology of the neoplastic epithelium. We aimed to compare genetic alterations in low-grade and high-grade regions of the same IPMN in order to identify molecular alterations underlying neoplastic progression. DESIGN: We performed multiregion whole exome sequencing on tissue samples from 17 IPMNs with both low-grade and high-grade dysplasia (76 IPMN regions, including 49 from low-grade dysplasia and 27 from high-grade dysplasia). We reconstructed the phylogeny for each case, and we assessed mutations in a novel driver gene in an independent cohort of 63 IPMN cyst fluid samples. RESULTS: Our multiregion whole exome sequencing identified KLF4, a previously unreported genetic driver of IPMN tumorigenesis, with hotspot mutations in one of two codons identified in >50% of the analyzed IPMNs. Mutations in KLF4 were significantly more prevalent in low-grade regions in our sequenced cases. Phylogenetic analyses of whole exome sequencing data demonstrated diverse patterns of IPMN initiation and progression. Hotspot mutations in KLF4 were also identified in an independent cohort of IPMN cyst fluid samples, again with a significantly higher prevalence in low-grade IPMNs. CONCLUSION: Hotspot mutations in KLF4 occur at high prevalence in IPMNs. Unique among pancreatic driver genes, KLF4 mutations are enriched in low-grade IPMNs. These data highlight distinct molecular features of low-grade and high-grade dysplasia and suggest diverse pathways to high-grade dysplasia via the IPMN pathway.

Genomic characterization of malignant progression in neoplastic pancreatic cysts.

Intraductal papillary mucinous neoplasms (IPMNs) and mucinous cystic neoplasms (MCNs) are non-invasive neoplasms that are often observed in association with invasive pancreatic cancers, but their origins and evolutionary relationships are poorly understood. In this study, we analyze 148 samples from IPMNs, MCNs, and small associated invasive carcinomas from 18 patients using whole exome or targeted sequencing. Using evolutionary analyses, we establish that both IPMNs and MCNs are direct precursors to pancreatic cancer. Mutations in SMAD4 and TGFBR2 are frequently restricted to invasive carcinoma, while RNF43 alterations are largely in non-invasive lesions. Genomic analyses suggest an average window of over three years between the development of high-grade dysplasia and pancreatic cancer. Taken together, these data establish non-invasive IPMNs and MCNs as origins of invasive pancreatic cancer, identifying potential drivers of invasion, highlighting the complex clonal dynamics prior to malignant transformation, and providing opportunities for early detection and intervention.

Intraductal Papillary Mucinous Neoplasms Arise From Multiple Independent Clones, Each With Distinct Mutations.

BACKGROUND & AIMS: Intraductal papillary mucinous neoplasms (IPMNs) are lesions that can progress to invasive pancreatic cancer and constitute an important system for studies of pancreatic tumorigenesis. We performed comprehensive genomic analyses of entire IPMNs to determine the diversity of somatic mutations in genes that promote tumorigenesis. METHODS: We microdissected neoplastic tissues from 6-24 regions each of 20 resected IPMNs, resulting in 227 neoplastic samples that were analyzed by capture-based targeted sequencing. Somatic mutations in genes associated with pancreatic tumorigenesis were assessed across entire IPMN lesions, and the resulting data were supported by evolutionary modeling, whole-exome sequencing, and in situ detection of mutations. RESULTS: We found a high prevalence of heterogeneity among mutations in IPMNs. Heterogeneity in mutations in KRAS and GNAS was significantly more prevalent in IPMNs with low-grade dysplasia than in IPMNs with high-grade dysplasia (P < .02). Whole-exome sequencing confirmed that IPMNs contained multiple independent clones, each with distinct mutations, as originally indicated by targeted sequencing and evolutionary modeling. We also found evidence for convergent evolution of mutations in RNF43 and TP53, which are acquired during later stages of tumorigenesis. CONCLUSIONS: In an analysis of the heterogeneity of mutations throughout IPMNs, we found that early-stage IPMNs contain multiple independent clones, each with distinct mutations, indicating their polyclonal origin. These findings challenge the model in which pancreatic neoplasms arise from a single clone. Increasing our understanding of the mechanisms of IPMN polyclonality could lead to strategies to identify patients at increased risk for pancreatic cancer.

Single-cell sequencing defines genetic heterogeneity in pancreatic cancer precursor lesions.

Intraductal papillary mucinous neoplasms (IPMNs) are precursors to pancreatic cancer; however, little is known about genetic heterogeneity in these lesions. The objective of this study was to characterize genetic heterogeneity in IPMNs at the single-cell level. We isolated single cells from fresh tissue from ten IPMNs, followed by whole genome amplification and targeted next-generation sequencing of pancreatic driver genes. We then determined single-cell genotypes using a novel multi-sample mutation calling algorithm. Our analyses revealed that different mutations in the same driver gene frequently occur in the same IPMN. Two IPMNs had multiple mutations in the initiating driver gene KRAS that occurred in unique tumor clones, suggesting the possibility of polyclonal origin or an unidentified initiating event preceding this critical mutation. Multiple mutations in later-occurring driver genes were also common and were frequently localized to unique tumor clones, raising the possibility of convergent evolution of these genetic events in pancreatic tumorigenesis. Single-cell sequencing of IPMNs demonstrated genetic heterogeneity with respect to early and late occurring driver gene mutations, suggesting a more complex pattern of tumor evolution than previously appreciated in these lesions. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Network Analysis of Protein Adaptation: Modeling the Functional Impact of Multiple Mutations.

The evolution of new biochemical activities frequently involves complex dependencies between mutations and rapid evolutionary radiation. Mutation co-occurrence and covariation have previously been used to identify compensating mutations that are the result of physical contacts and preserve protein function and fold. Here, we model pairwise functional dependencies and higher order interactions that enable evolution of new protein functions. We use a network model to find complex dependencies between mutations resulting from evolutionary trade-offs and pleiotropic effects. We present a method to construct these networks and to identify functionally interacting mutations in both extant and reconstructed ancestral sequences (Network Analysis of Protein Adaptation). The time ordering of mutations can be incorporated into the networks through phylogenetic reconstruction. We apply NAPA to three distantly homologous ß-lactamase protein clusters (TEM, CTX-M-3, and OXA-51), each of which has experienced recent evolutionary radiation under substantially different selective pressures. By analyzing the network properties of each protein cluster, we identify key adaptive mutations, positive pairwise interactions, different adaptive solutions to the same selective pressure, and complex evolutionary trajectories likely to increase protein fitness. We also present evidence that incorporating information from phylogenetic reconstruction and ancestral sequence inference can reduce the number of spurious links in the network, whereas preserving overall network community structure. The analysis does not require structural or biochemical data. In contrast to function-preserving mutation dependencies, which are frequently from structural contacts, gain-of-function mutation dependencies are most commonly between residues distal in protein structure.

Matching phenotypes to whole genomes: Lessons learned from four iterations of the personal genome project community challenges.

The advent of next-generation sequencing has dramatically decreased the cost for whole-genome sequencing and increased the viability for its application in research and clinical care. The Personal Genome Project (PGP) provides unrestricted access to genomes of individuals and their associated phenotypes. This resource enabled the Critical Assessment of Genome Interpretation (CAGI) to create a community challenge to assess the bioinformatics community's ability to predict traits from whole genomes. In the CAGI PGP challenge, researchers were asked to predict whether an individual had a particular trait or profile based on their whole genome. Several approaches were used to assess submissions, including ROC AUC (area under receiver operating characteristic curve), probability rankings, the number of correct predictions, and statistical significance simulations. Overall, we found that prediction of individual traits is difficult, relying on a strong knowledge of trait frequency within the general population, whereas matching genomes to trait profiles relies heavily upon a small number of common traits including ancestry, blood type, and eye color. When a rare genetic disorder is present, profiles can be matched when one or more pathogenic variants are identified. Prediction accuracy has improved substantially over the last 6 years due to improved methodology and a better understanding of features.

Evolution of Neoantigen Landscape during Immune Checkpoint Blockade in Non-Small Cell Lung Cancer.

Immune checkpoint inhibitors have shown significant therapeutic responses against tumors containing increased mutation-associated neoantigen load. We have examined the evolving landscape of tumor neoantigens during the emergence of acquired resistance in patients with non-small cell lung cancer after initial response to immune checkpoint blockade with anti-PD-1 or anti-PD-1/anti-CTLA-4 antibodies. Analyses of matched pretreatment and resistant tumors identified genomic changes resulting in loss of 7 to 18 putative mutation-associated neoantigens in resistant clones. Peptides generated from the eliminated neoantigens elicited clonal T-cell expansion in autologous T-cell cultures, suggesting that they generated functional immune responses. Neoantigen loss occurred through elimination of tumor subclones or through deletion of chromosomal regions containing truncal alterations, and was associated with changes in T-cell receptor clonality. These analyses provide insight into the dynamics of mutational landscapes during immune checkpoint blockade and have implications for the development of immune therapies that target tumor neoantigens.Significance: Acquired resistance to immune checkpoint therapy is being recognized more commonly. This work demonstrates for the first time that acquired resistance to immune checkpoint blockade can arise in association with the evolving landscape of mutations, some of which encode tumor neoantigens recognizable by T cells. These observations imply that widening the breadth of neoantigen reactivity may mitigate the development of acquired resistance. Cancer Discov; 7(3); 264-76. ©2017 AACR.See related commentary by Yang, p. 250This article is highlighted in the In This Issue feature, p. 235.

SubClonal Hierarchy Inference from Somatic Mutations: Automatic Reconstruction of Cancer Evolutionary Trees from Multi-region Next Generation Sequencing.

Recent improvements in next-generation sequencing of tumor samples and the ability to identify somatic mutations at low allelic fractions have opened the way for new approaches to model the evolution of individual cancers. The power and utility of these models is increased when tumor samples from multiple sites are sequenced. Temporal ordering of the samples may provide insight into the etiology of both primary and metastatic lesions and rationalizations for tumor recurrence and therapeutic failures. Additional insights may be provided by temporal ordering of evolving subclones--cellular subpopulations with unique mutational profiles. Current methods for subclone hierarchy inference tightly couple the problem of temporal ordering with that of estimating the fraction of cancer cells harboring each mutation. We present a new framework that includes a rigorous statistical hypothesis test and a collection of tools that make it possible to decouple these problems, which we believe will enable substantial progress in the field of subclone hierarchy inference. The methods presented here can be flexibly combined with methods developed by others addressing either of these problems. We provide tools to interpret hypothesis test results, which inform phylogenetic tree construction, and we introduce the first genetic algorithm designed for this purpose. The utility of our framework is systematically demonstrated in simulations. For most tested combinations of tumor purity, sequencing coverage, and tree complexity, good power (≥ 0.8) can be achieved and Type 1 error is well controlled when at least three tumor samples are available from a patient. Using data from three published multi-region tumor sequencing studies of (murine) small cell lung cancer, acute myeloid leukemia, and chronic lymphocytic leukemia, in which the authors reconstructed subclonal phylogenetic trees by manual expert curation, we show how different configurations of our tools can identify either a single tree in agreement with the authors, or a small set of trees, which include the authors' preferred tree. Our results have implications for improved modeling of tumor evolution and the importance of multi-region tumor sequencing.

A probabilistic model to predict clinical phenotypic traits from genome sequencing.

Genetic screening is becoming possible on an unprecedented scale. However, its utility remains controversial. Although most variant genotypes cannot be easily interpreted, many individuals nevertheless attempt to interpret their genetic information. Initiatives such as the Personal Genome Project (PGP) and Illumina's Understand Your Genome are sequencing thousands of adults, collecting phenotypic information and developing computational pipelines to identify the most important variant genotypes harbored by each individual. These pipelines consider database and allele frequency annotations and bioinformatics classifications. We propose that the next step will be to integrate these different sources of information to estimate the probability that a given individual has specific phenotypes of clinical interest. To this end, we have designed a Bayesian probabilistic model to predict the probability of dichotomous phenotypes. When applied to a cohort from PGP, predictions of Gilbert syndrome, Graves' disease, non-Hodgkin lymphoma, and various blood groups were accurate, as individuals manifesting the phenotype in question exhibited the highest, or among the highest, predicted probabilities. Thirty-eight PGP phenotypes (26%) were predicted with area-under-the-ROC curve (AUC)>0.7, and 23 (15.8%) of these were statistically significant, based on permutation tests. Moreover, in a Critical Assessment of Genome Interpretation (CAGI) blinded prediction experiment, the models were used to match 77 PGP genomes to phenotypic profiles, generating the most accurate prediction of 16 submissions, according to an independent assessor. Although the models are currently insufficiently accurate for diagnostic utility, we expect their performance to improve with growth of publicly available genomics data and model refinement by domain experts.

Exome sequencing identifies frequent inactivating mutations in BAP1, ARID1A and PBRM1 in intrahepatic cholangiocarcinomas.

Through exomic sequencing of 32 intrahepatic cholangiocarcinomas, we discovered frequent inactivating mutations in multiple chromatin-remodeling genes (including BAP1, ARID1A and PBRM1), and mutation in one of these genes occurred in almost half of the carcinomas sequenced. We also identified frequent mutations at previously reported hotspots in the IDH1 and IDH2 genes encoding metabolic enzymes in intrahepatic cholangiocarcinomas. In contrast, TP53 was the most frequently altered gene in a series of nine gallbladder carcinomas. These discoveries highlight the key role of dysregulated chromatin remodeling in intrahepatic cholangiocarcinomas.

Network models of TEM ß-lactamase mutations coevolving under antibiotic selection show modular structure and anticipate evolutionary trajectories.

Understanding how novel functions evolve (genetic adaptation) is a critical goal of evolutionary biology. Among asexual organisms, genetic adaptation involves multiple mutations that frequently interact in a non-linear fashion (epistasis). Non-linear interactions pose a formidable challenge for the computational prediction of mutation effects. Here we use the recent evolution of ß-lactamase under antibiotic selection as a model for genetic adaptation. We build a network of coevolving residues (possible functional interactions), in which nodes are mutant residue positions and links represent two positions found mutated together in the same sequence. Most often these pairs occur in the setting of more complex mutants. Focusing on extended-spectrum resistant sequences, we use network-theoretical tools to identify triple mutant trajectories of likely special significance for adaptation. We extrapolate evolutionary paths (n = 3) that increase resistance and that are longer than the units used to build the network (n = 2). These paths consist of a limited number of residue positions and are enriched for known triple mutant combinations that increase cefotaxime resistance. We find that the pairs of residues used to build the network frequently decrease resistance compared to their corresponding singlets. This is a surprising result, given that their coevolution suggests a selective advantage. Thus, ß-lactamase adaptation is highly epistatic. Our method can identify triplets that increase resistance despite the underlying rugged fitness landscape and has the unique ability to make predictions by placing each mutant residue position in its functional context. Our approach requires only sequence information, sufficient genetic diversity, and discrete selective pressures. Thus, it can be used to analyze recent evolutionary events, where coevolution analysis methods that use phylogeny or statistical coupling are not possible. Improving our ability to assess evolutionary trajectories will help predict the evolution of clinically relevant genes and aid in protein design.